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An algebraic approach to the neutrino propagation in dense media is presented. The Hamiltonian 
describing a gas of neutrinos interacting with each other and with background fermions is written in 
terms of the appropriate SU(N) operators, where N is the number of neutrino flavors. The evolution 
of the resulting many-body problem is formulated as a coherent-state path integral. Some commonly 
used approximations are shown to represent the saddle-point solution of the path integral for the 
full many-body system. 
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I. INTRODUCTION 



O: 

' ■ Neutrino propagation in dense matter is encountered near the core of a core-collapse supernovae [_ ^ , 

in the Early Universe 0, [Tol |. and possibly in the gamma-ray bursts In particular, neutrino interactions play 

^ . a crucial role in core-collapse supernovae. Neutrino oscillations in a core-collapse supernova differ from the matter- 
■ enhanced neutrino oscillations in the Sun as in the former there are additional effects coming from both neutrino- 
CN . neutrino scattering 0, 0, ^| and antineutrino flavor transformations 0| . Exact integration of the neutrino evolution 
' equations with such terms in the supernova environment turns out to be a very difficult problem. A "mean-field" type 
approximation was proposed in Ref. 01 which was adopted in exploratory calculations of the conditions for r-process 
nucleosynthesis in supernovae 0] . Analytical solutions were investigated in the limitin g ca ses where off-diagonal 
\ terms dominate . In addition, various collective effects were explored 0, 0, 0, 0, ITMTl^ l2Cj | . 
■~>«^ I The purpose of this article is to formulate the problem of neutrino propagation in dense media algebraically. An 
' algebraic formulation of the problem should make hidden symmetries evident and also provide a framework to look 
for exact solutions or systematic approximations. In the next section we introduce the second-quantized formalism 
Q ' and the appropriate SU(2) algebra for two flavors. In this section we also show that the evolution operator for the 



standard MSW problem of neutrinos, mixing with each other and interacting with background electrons, can be 



^ written down exactly using an algebraic ansatz. In this and subsequent sections we utilize a path-integral approach to 
the underlying many-body problem, the details of which are sketched in the Appendix. In Section III, we introduce 
^ the neutrino-neutrino interaction and elucidate its algebraic nature. These arguments are expanded into the situations 
where antincutrinos arc present in Section IV. The full problem with three flavors of both neutrinos and antincutrinos 
is discussed in Section V. Finally, a discussion of our results in Section VI concludes the paper. 

II. SECOND-QUANTIZED FORMALISM 



In our presentation we find it convenient to use the second-quantized form [2l|. For simplicity in this section we 
assume that only two flavors of neutrinos mix, which we take to be the electron neutrino, i^e, and a combination 
of muon and tau neutrinos, which we denote by v^- (In the limit ^13 = 0, one particular combination of the mu 
and tau flavors decouple and the results in this section become exact with i^x being the combination orthogonal to 
the decoupled one |23|). We first consider a situation where there are no antincutrinos. (We relax both of these 
assumptions in the subsequent sections). The Hamiltonian describing mixing and interaction of the neutrinos with 
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the background electrons is given by 

^"^IJ "^'^ ^"""^^ " " cilip)aeip)) 

+ IJ'^^'P^ sm2e{al{p)aeip) + al{p)a4p)), (1) 

where a\{p) and aj.(p) are the creation operators for the left-handed i^e and i^x with momentum p, respectively, and 
ae{p) and a^ip) are the corresponding annihilation operators. In Eq. (QJ, ^/2GFNe{x) is the Wolfenstein 

potential describing the interaction of neutrinos with electrons in neutral, unpolarized matter 23]; = n^- — ne+ is 
the net electron density; is the vacuum mixing angle; and 5m^ = m\ — m\. In writing Eq. (^) we omitted a term 
proportional to the identity (this includes the other Wolfenstein potential Vn{x) = — (1/\/2)G_f-/V„(x) describing the 
neutral-current interaction of neutrinos with neutrons). The presence of the Mikheev, Smirnov, Wolfenstein (MSW) 
resonance [23,13113 manifest in the first term. 

The inherent SU(2) symmetry of the problem can be implemented by the operators 

J+{p)^al{p)ae{p), J-{p) ^ al{p)ax{p), Jo{p) ^ ^ {al{p)ax{p) - al{p)ae{p)) , (2) 
which satisfy the commutation relations 

[J+(p), J-iq)] = 25\p - q)Joip), [Jo(p), J±iq)] = ±SHp - q)J±{p). (3) 

These equations describe as many commuting SU(2) algebras as the number of distinct values of the neutrino momenta 
p permitted by the appropriate physical situation. Each J{p) is realized in the j — 1/2 representation due to the 
fermionic nature of neutrinos. The global SU(2) operators, 

j± = y" dv±(p), Jo = y" dVob), (4) 

also satisfy the SU(2) commutation relations and play an important role in the problem. Representation of the algebra 
spanned by the operators J of Eq. is obtained by adding N different copies of SU(2) each with j = \/2 where N 
is the number of allowed values of neutrino momenta. 

In terms of the SU(2) generators of Eq.Q, we can write the neutrino Hamiltonian in Eq. ^ as follows: 



COS 26*70 (p) + ^ sin 26* {J+{p) + J_(p)) 



\/2Gf j d^pN^Mp). (5) 



Note that in the last term of Eq. jsj we kept the electron density inside the momentum integral since it depends 
on the direction of the neutrino momentum (neutrino traveling in different directions will, in principle, see different 
electron densities), but of course not on the absolute value. 

In general we are interested in finding the operator [/, describing the evolution of the system: 

dU 

^^ = H.U. (6) 

In most cases, however, the physical interest is in finding the evolution of a particular state of the system. One may 
start, for example, with a state in which all permitted electron neutrino states are occupied and all states are 
empty i.e., 

i'^) = n«^(p)io)' 

where |0) is the particle vacuum. V denotes the set of all allowed neutrino momenta. The state \4>) is annihilated by 
the operators J-{p) for all p and also by JL. It is the lowest-weight eigenstate of the representation with j = N/2 of 
the global SU(2) operators Jq and J±. 

The evolution operator can be found by employing the unitary ansatz 

U^expf [ d^pT{p,t)J+{p)) exp ( [ c;3plog(l + |r(p, i)!') Jo(p)) exp (- [ d^pr* {p,t)J-{p)] . (8) 
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Here T{p,t) is a function to be determined by substituting this ansatz into Eq. JBJ. One can differentiate U of 
Eq. using the operator chain rule. Differentiation introduces the operator Jo{p') between the first and second 
exponentials, and the operator J-{p') between the second and third exponentials. Since, for example, Jo{p') does not 
commute with the first exponential, to write this operator before the first exponential, one needs to introduce the 
identity operator in the form exp (J d'^pT{p,t)J+{p)) exp (— J (PpT{p,t)Jj^.{p)) after the operator Jo{p') and use the 
identity 

exp(0i)02 exp(-Oi) = 0-2 + [Oi, O2] + l^Pi, [d, O2]] + • • • , 
valid for any two arbitrary operators Oi and Oi. Moving all such terms to the left of the exponentials we find 



dU 



,3 ( ■ Tip,t)T*{p, t) ~ t{p, t)f*{p, t) 



u, (9) 

where the dot denotes derivative with respect to time. Eq. (j^J is satisfied if r(p, t) obeys the equations: 



T(p,t)T*(p,t) - T(p,t)f*(p,t) Sm^ r- , , 



l + |T(p,t)|2 2p 

and 

T{p,t) , T*{p,t) 1 Sm 



sin 26*, (11) 



l + |T(p,t)|2 l + |r(p,i)|2 2 2p 

i.e. equations for different neutrino momenta decouple as expected. In order to satisfy the condition U{t = 0) = 1, 
t{p, t = 0) should be zero for all p. Using Eqs. H1U|) and (|ll|l . one can show that t{p, t) is the solution of the following 
nonlinear (Riccati-type) first-order differential equation: 

ifip,t) = ^^-^sm2e{l~T{p,t)')+(^^cos29~V2GFN,)j T{p,t). (12) 

r(p, t) can be interpreted as the ratio of the one-body neutrino wave functions 

np, t) = I ( ,N ■ (13) 

To see this, one only needs to substitute Eq. (|13|) into Eq. (|12f) together with the normalization condition 

IV'eP + I^.P^l. (14) 

This way, it becomes clear that Eq. (|12|l is equivalent to the Schrodinger equation 



._9 / Ve A 1 Acos26i A sin 20 \ ( '^^'^ 

^dt\i^x)^2\ Asin26' -A + A cos 26* i i 



(15) 



where we defined 



A = A = V2GfN,. (16) 

2p 

For T{t = 0) to be zero, we need the initial conditions ^|Jx{t — 0) — and ipe{t = 0) = 1. Eq. H15(l is nothing but 
the standard MSW evolution equation for two flavors. Note that our ansatz for the evolution operator provides the 
"coset-space" formulation of this problem (cf. Eq. (|12|l with the formulation in the Appendix of Ref. ,2^]). The 
initial state given in Eq. 0, then evolves into the state 



exp (~IJ^ d'plog(l + \t{p, t)\^)^ exp (^J^ d^pT{p, t)J+ 



ip) 10), (17) 



where T{p,t) is obtained by solving Eq. (|12|) . Both the initial state in Eq. (TJ and evolved state in Eq. (|17|l are 
normalized to unity. The evolved state in Eq. 1171) is the standard normalized SU(2) coherent state. 
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III. NEUTRINO-NEUTRINO INTERACTIONS 



In some astrophysical environments, such as supernovae and the Early Universe, where neutrino density can become 
very large [23|, the Hamiltonian in Eq. is no longer sufficient. Neutrino self interactions must also be taken into 
account. In most cases one is mainly interested in the forward scattering of neutrinos from other neutrinos. These 
are described by the Hamiltonian 



d^pd^qRpq [al{p)ae{p)al{q)ae{q) + a^^{p)ax{p)a\{q)a^{q) 



V2V . 

al{p)ae{p)al{q)a^{q) + al{p)ax{p)al{q)ae{q)] (18) 



where we defined 



= (l-cos?9p,). (19) 



Here, "dpq is the angle between the momentum directions of the neutrinos with momenta p and q and we used 
box quantization conditions for a box with volume V. In forward scattering, neutrinos exchange their momenta by 
exchanging a Z boson with momentum p — q or a Z boson with zero momentum. The coefficient 1 — cos -dpg ensures 
that the neutrinos which move parallel to each other do not undergo forward scattering. 

Numerical studies indicate that the addition of a neutrino background results in very interestin g ph ysical effects 
such as coherent flavor transformation. For details see, for example, Refs. H S H H and [H IHlll, IH 13 • The 
major difficulty in studying the effects of a neutrino background on flavor evolution is the inherent nonlinearity of the 
problem. On the other hand, the contribution from H^,^ is relatively easy to study in the algebraic formalism. First, 
one needs to write H^^, given in Eq. H18|) in terms of the SU(2) generators of Eq. Omitting a term proportional 
to the identity, we find that 



H.. = / d'p d'q Rpg J(p) • J(g). (20) 



Here, 3{p) ■ 3{q) is defined as 



3{p) ■ 3{q) = Jo{p)Mq) + \j+{p)J-{q) + \j-{p)J+{q). (21) 

This expression is very similar, but not identical, to the global SU(2) Casimir operator. In fact, Hyi, commutes with 
the global SU(2) operators of Eq. Q 

[i?.., J.] = (22) 

although the Hamiltonian does not: 

[i?.,Ji]^0. (23) 

On the other hand, the last term of Hi, in Eq. is proportional to Jq for a constant electron density. Thus if 
5m? were zero the evolution problem for + H,ji, could have been solved in the j7-basis for a constant electron 
density. We see that, although H^i, presents a difficulty in numerical studies, it is relatively easier to work with in 
this algebraic approach. The noncommutativity of the first integral in Eq. (|SJ) with generators of the total angular 
momentum operators J , however, is the salient difficulty of the algebraic approach to the problem. 

One possible approach to the problem of finding the evolution operator with Hi, + H^i, is to seek a path integral 
representation. One can appropriately use the SU(2) coherent states 

|z(t)) =iVcxp (^l^d^pz{p,t)J+{p)ycj,). (24) 

Here |0) is the state defined in Eq. (O and the normalization constant N is given by 

7V = exp(^-i^d3plog(l + |z(p,i)n) . (25) 
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Clearly when z{p,t) = T{p,t) the state in Eq. H24() becomes the exact solution of the evolution with alone. Path 
integral representation of the matrix element of the evolution operator calculated with H^, + H^,^ between two states 
\z{ti)) and \z'{tf)) is given by 

{z'{tf)\U\z{U))^ I i?[z,z*]e^^[^'^*I (26) 

where the path integral measure is 

D[z,z*] = hm e-^i:lJ.'^^'°H(i+WP^*^)0 TT TT ^l^I^. (27) 

a=l peV 

Here, the exponential factor arises because the SU{2) coherent states are overcomplete and therefore require a weight 
function in the resolution of identity. A detailed derivation can be found in the Appendix. In the above formula dz"^ 
refers to dzdz* . The action functional is given by^ 

S[z,z*]^ j^' dt{i^^-H,-H,,)-i\og{z'{tf)\z{tf)). (28) 

The leading contribution to this path integral comes from the stationary path \z{t)) which minimizes the action 
functional. As in any variational problem, stationary path can be found by solving the Euler-Lagrange equations 

where 



L[z,z*]^{i^-H,-H,,) (30) 
ot 

which plays the role of the Lagrangian. To solve the Euler-Lagrange equations, however, we first need to calculate 
the Lagrangian as a function of z and z* . This can be done by substituting from Eq. ISJ and H^i, from Eq. H2()|l 
into the Lagrangian of Eq. H3U|) and using the expectation values 

{j+{p)) ^ {j^{p)Y = "]^^:\, , (Mp)) = !7|"|^'!i'' (31) 

1 + \z{p,t)\^ 2 1 + \z{p,t)\'' 

which are valid for p E V. When p ^ V, however, these expectation values are equal to zero. When p ^ q coherent 
states also obey 

{Jaip)Jliq))^{Ja{p)){Jliq)) (32) 
for a — 0,zt. Using above formulas and Eqs. |(SJ), (|30|l and H2U|) together with 

'dt^ V ^l + k(p)| = 



i^^J-^ldpi^^, (33) 



we obtain the following expression for the Lagrangian: 



3 Bip){l - |z(p)|2) + 2B,M<P) + 2B,e{p)z*ip) 

l^[Z,^ J = - -VZUr_F J " 

Here we defined 



/ ,o z(p,t)z*(p,t) 1 f f6m^ „^ /-^ , \ 1 - |z(p,t)|2 
U ^ i d^p ,\ + - / d^p — — cos26i-V2Gi.A^e' 



V 



l + |z(p,i)P ■ 2U ^\2p ^ V l + k(p,t)P 



l/d3p^sin2^-(^'^)+/*t'^ 
2Jv 2p l + |z(p,i)P 



d>— — sin26> ' — (35) 



1 We use the shorthand notation (O) = {z{t)\0\z{t)) 
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-B(P) = / 0^ g-^pg i , -^^eb) = Sg^r = 2—— / d qRpg (36) 



and 

^W- p. '^"'"'1 + 1^(9, OP' ^:.ev., - ^e. - - p. " ^^^^ 1 + |2 ' 

Euler-Lagrange equations which follow from this Lagrangian are 

ii{p,t) = P{p,t) - a{p,t)z{p,t) ~ p*{p,t)z{p,tf (37) 

and its complex conjugate. The coefficients a and (3 are given by: 

StjIj^ 1 StjIi^ 1 

= -^cos20 + V2GfiVe + S(p), ^(p, t) = sin 20 + -B,e(p). (38) 

Zp 2 2p 2 

As in the previous section, ^(p, i) can be interpreted as the ratio of the one-body neutrino wave functions V'e, V'a;- 

= (39) 

Ve{P,t) 

with the normalization condition 

|^eP + |V'xP = l. (40) 

If we substitute Eq. ^ into Ea. (|77|) . we see that Eq. (|37|l is equivalent to the Schrodinger equation 



._a / -0e \ 1 / A + B- Acos26' Sea; + A sin 26* \ f ipe 
^dt\^Jx)~2\ B^e + Asin26i -A - S + A cos 26* j 

where A and A are defined in Eq. (|16|l . B and Be^; are obtained by substituting Eq. I|39|l into Eq. I|36|) : 



(41) 



^ = /_ [l^e(9,t)|' - 1^.(9, i)n , (42) 



Be. = ^ V'e (g, t)r, iq,t)- (43) 

Eg. H41() is commonly used to study neutrino propagation with neutrino- neutrino interactions (see Refs. Q through 
PJj). Here we illustrated that it is not an exact result, but represents the saddle point solution to the many-body 
problem. 

There is an alternative method to obtain the approximate results outlined above. This method, typically employed 
in the random-phase approximation approach to many-body problems, approximates the product of two commuting 
arbitrary operators Oi and O2 as 

0102 ^ Oi(^io2io + {^\oi\oo2 - {^mom2\o, (44) 

provided that the condition 

(ci6ia2io = (ciOiio(eio2io (45) 

is satisfied. In Eq. (|44|l the state |^) must be appropriately chosen so that Eq. H45|l is satisfied. Since the quantity 
Rpq is zero when the operators do not commute (i.e., when p = q), one can apply this approximation technique to 
the quadratic element 3{p) ■ 3{q) which appear in the neutrino- neutrino forward scattering Hamiltonian H^i, given in 
Eq. H2UI) . It follows from Eq. (|32|l that the SU(2) coherent states satisfy the condition stated in Eq. H45I) . Using the 
symmetry Rpq = Rqp, one can then replace H^^, by 



H.. ~ f <fpd\Rpq (joip){Mq)) + ^J+(p)(J-(<z)) + ij_(p)(J+(<z)) ) , (46) 



where the expectation values are calculated using the SU(2) coherent states (the variable z{t) of these coherent states 
must be obtained in a self-consistent way). In writing Eq. H46|) we omitted a term proportional to identity. The 
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total Hamiltonian for the two neutrino flavors propagating in the presence of electrons and other neutrinos in the 
background becomes approximately 



+n,, ^ J d^p {^-^ cos 26 - V2GFNe - B{p)^ Mp) 

+ J d^P Qsin2(?+ is.eb)^ J+{P) + J d^P Q sin + ^^^.(p)^ J 



(47) 



-ip)- 



Here B{p), Bex{p) and B^dp) are given in Eq. and we used the expectation values given in Eq. H31|) . Since 

the above approximate Hamiltonian is linear in the SU(2) generators we can use the same ansatz for the evolution 
operator of the system as we did in section H, i.e. 



d^pz*{p,t)J-ip) 



(48) 



(49) 



U = exp (^J d''pz{p,t)J+{p)^ exp (^J d3plog(l + \z{p,t)\^)Jo{p)^ exp J 
It is straightforward to show that, substituting this ansatz into the equation 

dt 

yields nothing but Eq. (|37|) for z{p, t). It is also straightforward to generalize this linearization scheme to the situations 
where antineutrinos and more flavors are present. 

Corrections to the path integral will naturally arise as a result of the deviations from the classical (i.e., stationary) 
path given by Eq. H37I) or Eq. (|41|l . To calculate the effect of small deviations one can carry out a Taylor expansion 
of the action around the classical path. It is helpful to introduce the following nonlinear transformation in order to 
eliminate the exponential factor in the path integral measure in Eq. (|27(l : 



Pip.t) 



z{p, t) 



V^ + \zip,tW 



P*{p,t)^ 



z*{p,t) 



^1 + \z[pM 



(50) 



The Jacobian resulting from this change of variables conveniently cancels out the exponential factor. It is straight- 
forward to show that 



N 



The path integral then takes the form 



d'^z{p, tg) 

2m 



N 



nn2! 



d^l3{p,tg) 

2TTi 



{z'{tf)\U\z{U)) 



N 



dpjp, ta)d(3*{p, tg) ^^s\p,(}* 

N^oo ITT 
a=l p<£-p 



JL- n n 



(51) 



(52) 



where the action functional, written in terms of the new variables is as follows: 



S[f3,f3* 



dt 
1 6m'^ 



dp 



P{p,t)P*ip,t)~P(p,t)p*ip,t) , 1 fSm 



2i 



2p 



cos29-V2GfnA {l~2\l3{p,t)\ 



2 2p 

1 V2Gf 

2 V 



sin 26* 



(^p{p,t)+p*{p,t)^l-\/3{p,tW 
d^pd\Rp, [(|/3(p,t)P + |/3(<z,0l' - 2|/3(p,i)|2|/3(q,t)|2) 



(53) 



+ V(l - mP.W) (1 - {P{p,t)P*{q,t)+P*{p,t)P{q,t)) 
One obtains the following equation of motion (and its complex conjugate) from the variation of this action: 



iP{p,t) 



2p 
6m? 



coa2e - V2GfN, 
V2Gf 



V2Gi 
V 



2p 



■ sin 20 



2p 



■ sin 29 - 



V 

V2Gf 
V 



d^qRpW^-\PiQ,tWPiQ,t) 



d^qRp,y/T^^P{^(3{q,t) 



d^qRpg{l-2\f3iq,t)\^)j(3{p,t) 

' mp,tr~2 

2Vl-|/3(p,t)P 
P{p,tf 



(54) 
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This classical path is the same as the one given in Eq. H37(l . This can be shown directly by substituting the 
transformation in Eq. H5U|1 into Eq. (j^TJ. We will denote the classical path as f3ci- Since the first order variations are 
zero on the classical path by definition, the Taylor expansion around Pel yields 

S[(3,(3*] ^ S[f3ci,P*ci] + Uf^~M^ (y^) iP-M (55) 



2 r-c., ygpgp 



Here {■ ■ ■)ci indicates that the derivatives are to be calculated on the classical path Pel and 
{ww) ~ ^'^'^ ^ hand notation for the matrix product 

E E (P(P^ - MP, t.)f - MP, t^)) , (56) 

and similarly for the other terms. The sums over p and q run through the allowed momentum modes and the sums 
over k and m run from 1 to which is the number of time intervals we introduced in path integral. N ^ oo limit 
should be taken as explained in the Appendix. For small deviations from the classical path, one can ignore the higher 
order terms in the expansion (|55|l and substitute it in Eq. 152|1 : 

(57) 

where we defined P = P — Pel- The classical action S[Pci, P*i] is taken out of the integration since it does not depend 
on p. The lowest order quantum corrections are captured by the Gaussian integral in Eq. (|57|l . The result of the 
integration is 

{z'{tf)\U\z{U)) = lim (i7r)^+^^=^^^^^=. (58) 
Here P denote the number of allowed momentum modes. The matrices M and L are given as follows: 

^' = ^ (^.(,,J)l(,,t.)) , ''^^^ k,q^-^)-\ (,,(,,,,)X(,,t„)) , 

L{p,k,q,m) 2 \Sx{p,tk) Sy{q,t,n) ) ^^^^ 

where x = {P + P*)/2 and y = [P — P*)/2i. The fundamental difficulty involved in the calculation of the determinant 
in Eq. H58fl is the existence of non-diagonal terms in the matrices M and L. These terms are generated by the 
/ d^p (Pq Rpq . . . integral in the action given in Eq. H53|l . A complete analysis of these determinants is beyond the 
scope of the present paper and will be given elsewhere. 

IV. NEUTRINO-ANTINEUTRINO INTERACTIONS 

Environments such as core-collapse supernovae and the Early Universe contain copious amounts of antineutrinos 
as well as neutrinos. When antineutrinos are added to the picture the effective flavor evolution Hamiltonian becomes 

H^H^ + H^ + H^^ + H^^ + H^^. (61) 

Here, H,y and H^^ are given in Eqs. ^ and ((TH|) . iJp and iJpi> are the same as and H^^, respectively except that 
the neutrino operators a and are replaced by the antineutrino operators b and and the sign of the A^e term is 
reversed. The neutrino- antineutrino forward scattering Hamiltonian Hi,p is given as follows: 



V2G 



F 



d^pd^qRpq [alip)aeip)bl{q)beiq) + aUp)a4p)bl{q)b4q) (62) 



V 

aUp)ae{p)bUq)be{q) + al{p)a,ip)bl{q)b,iq)] . 
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In addition to the neutrino operators defined in Eq. (0), we now define the antineutrino SU(2) operators^ 

J+ip)^blip)b,ip), J_ip)^blip)KM, Mp)^l{blipMp)-bUp)bAp))- (63) 
These operators also obey SU(2) commutation relations 

[J+(p), J_(q)] = 2S{p-q)Mp), [Mp),J±{q)] = ±5{p~q)J±{p) (64) 
and they commute with the neutrino operators 

[Mp),Mq)]=0. (65) 
Written in terms of Ji{p) and Ji{p), the Hamiltonian in Eq. (|61|l takes the following form: 

H^+Hi,+ H^^ + iJpp + (66) 
= j d?P {^-^ cos2^? - x/2Gf iVe^ Up) + i y Sp^-^ sin 20 (J+(p) + J_(p)) 

d^P f - — cos20 - ^GfNA Jo{p) + -j d^P^ sin 20 (J+(p) + J_(p)) 

' 1 d^pdSqi?^, (J(p) . 3{q) + 3{p) ■ J(g)) 



V 
V2Gi 



d^pd^qRpg {-2Jo{p)Jo{q) + J+{p)J-{q) + J-{p)J+{q)) 



V 

where a term proportional to identity is omitted. 

As we did in the previous section, we introduce the coherent states'^ 

\z{t),z{t)) = NNe^v d'p^(p-.t)j^ip)J^ d^pHP,t)JAP)^^y (g7) 

Here \(f>) is analogous to the state defined in Eq. ((TJ, i.e., all permitted Ve and i^e states are occupied and all other 
neutrino flavor states are empty: 

10) = n «e(p) n ^e(p)IO). (68) 

pev per 

In the above formulas, V and V denote the set of all allowed momentum modes for neutrinos and antineutrinos, 
respectively. In what follows, we will drop the symbols V and V from the formulas by adopting the convention that 
the non-overlined quantities such as p, q will take values in V whereas the overlined quantities such as p, q take values 
in V. 

The constants 

N - exp (^-i J d^p\og{l + \z{p,t)\^)^ , TV - exp (^-^ J d^plog{l + |z(p,i)|')) (69) 
in Eq. 1)67(1 normalize the coherent states: 

(z,z|z,z)-l. (70) 
A path integral representation of the evolution operator can be given in terms of these coherent states as 

{z'{tf ),z'{tf )\ U \z{U),z{U)) = I D[z, z\ z, -z*YS[z..'r..r\ (71) 



^ Time-reversal invariance requires that the order of the flavors is flipped in the deflnition of the antineutrino algebra as compared with 

the definition of the neutrino algebra in Eq. J^J- 
^ In these formulas z is an independent complex number, not the complex conjugate of z. Complex conjugates will be denoted by a star 

such as z* and 2*. 



10 



where the action functional S[z, z* , z, z*] is'' 

S[z, z*,z, ^*] = J^' dt{i^^ - {H, +H^ + H,, + H,^ + H,^)) - i ln(z'(i/), z' {tf)\z{tf),z(tf)) . (72) 

Once more we can find the stationary path by solving the Euler-Lagrange equations for 

L[z, z\-z, -z*\ - - {U, + H-, + H,, + + H,,-,)). (73) 
at 

The linear terms in the Lagrangian can be calculated using 

which are valid for p and p CzV. When p oi p the expectation values are zero. To calculate the quadratic 
terms we use the identities 

{Ja{p)jl{<l)) = {Ja{p)){jl{q)). {Ja{p)jlm = ( Ja (p)) ( 4 (g)) , {JaipUUq)) - {Ja{p)) {jHq)) ■ (76) 

Here a — 0,± and we assumed p ^ q and p ^ q. Also note that 

^ dt' \J ^l+|z(p)|2 J ^l+|z(p)|2; 

Using Eqs. (|74() - (|77|l Lagrangian can be found as follows: 



1 /■ ,3 i?(p)(l ~ k(p)P) + 2BUp)z{p) + 2BUp)z*{p) 



1 f^3-BiP)a - \m\^) + ^B^eiPMp) + 2B,^{p)-z*{p) 



Here L^, is the Lagrangian given in Eq. H35I) and L^, is the same as L^, except that we substitute z[p) in place of z{p) 
and change the sign of N^.- In the above equation we also defined 



and 



( \ ( ^ O^l^L ( rfi R ^(g.^) n ^Gp /• ,3_p Z*{q,t) 



Equations of motion which result from this Lagrangian are as follows: 

izip, t) = I3(j}, t) - aip, t)z(j}, t) - p* {p, t)z^{p, t), (81) 

ii{p, t) = -P* ip, t) + a(p, t)z{p, t) + P{p, t)z\p, t) (82) 



* Here the expectation values are calculated using the states in Eq. 1671 : 

{O) = {z{t),z(t)\0\z(t),z(t)) 
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and the complex conjugates of these equations. The coefficients a and /3 are given by: 



a(p, t) 



2p 



cos 20 + V^GpNe + B{p) f3{p, t) 



2 2p 



sm29 + Be^{p). 



(83) 



The coefficients a and f3 are the same except that the sign of the Sm^ is changed. 

We can again write the parameters z and z in terms of the one-body neutrino and antineutrino wave functions ipe, 
ipx^ ipe and ij^x as follows: 



z{p, t) 



and z{p, t) 



Aip, t) 

We also have the normalization conditions 

\Af + \^xf = l and l^el' + IV^xP = 1. 
Substituting Eq. ^ into Eqs. ^ and ^ we see that they are equivalent to the Schrodinger equations 



.0 ( ^eip,t) \ _ I f A + B - Acos2e Bex + Asm2e 



and 



at V ^Ap,t) 



dt \ tpx{p,t) 



1 / A + B + A cos 261 Bex- A sin 29 
Br,. - A sin 261 -A- B- A cos 261 



tpeiP^t) 

Mp,t) J ' 



-06 (P, t) 

■4'x{p,t) 



In the above formula B and B^x are obtained by substituting Eq. (|85|) into Eqs. H79|) and (|80f) . i.e., 



B = 



V2Gi 
V 



V2Gi 
V 



d^qRpg{\Mq,t)\^-\MQ>tW 



Bex = B*e — 2 



V2Ch 
V 



d^qRpqiPe{q,t)^l^l{q,t)-2 



V2Gi 
V 



d^qRpgilJeiq,t)iJj*{q,t). 



(84) 



(85) 



(86) 



(87) 



(89) 



V. THREE NEUTRINO FLAVORS 

In this section we generalize our formalism to three neutrino flavors, i.e. when ^13 ^ 0. SU(3) symmetry of the 
neutrinos can be represented by the operators 



Tt] (p) = 4 {p)aj {p) and Ty- (p) = b] {p)b, {p) 



(90) 



where i and j run over the flavor indices e, /i and r. These operators generate orthogonal SU(3) algebras: 

[Tij{p),Tki{q)] ^ S{p - q) {SkjTuip) - SuTkjip)) , (91) 
[T^J{p),Tkl{q)] = -S{p~q) {6kjTu{p) - Sufkj{p)) , 
mj{p),fki{q)] = 0. 

The effective Hamiltonian describing the propagation of neutrinos in a background of electrons is given by p9l | 

(92) 



H = ^ d^p [(7jj (p) + Lj.,j (p)) T^j (p) + (p) - uJ^j (p)) T^j {p)] 
+ ^1 d'pd\RpgJ2{T.,{p)~T,,{p)){T,,{q)-f,,{q)). 



In this Hamiltonian, the coefficients ^ij{p) are the elements of the symmetric, traceless matrix F which is given below; 

\ 



r = hi]) = ^Q23Qi3Qi2 



/ -<5m2l-'^™3l 
2p 



V 



2p 








2p 



*9l2*9l3*923 



(93) 
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Here Q12, Q13 and Q23 are neutrino mixing matrices in vacuum: 

/I \ / Ci3 5l*3\ / C12 Si2 0\ 

Q23Q13Q12 - C23 ^23 1 -^12 C12 , (94) 

VO -^23 C23/ \-Si3 C13/ V 1/ 

where C13, etc. is the short-hand notation for cos0i3, etc. Since 5*13 may be muhipUed by a phase, we exphcitly 
indicated its complex conjugate. Note that individual matrices, not their matrix elements, are called Q23, Qi3, and 
Q12, respectively in Eq. 

The terms in the Hamiltonian which are proportional to ^ij represent vacuum oscillations of the neutrinos. The 
coefficients ujij are real and symmetric. They are the elements of the diagonal, traceless matrix 57 given below: 



214 

~Vc 
~Vc 



n = (uj,,) = i I ~Vc I . (95) 



Here Vc is the Wolfenstein potential described earlier. 

A path integral formula for the evolution operator of the system can be constructed using SU(3) coherent states 
which are given by 

These coherent states are defined with respect to the reference state {(p) which is defined as in Eq. H68() . The 
normalization constants N and TV are given by 

7V = expf-i / d'p\og{l + \z^ip,t)\^ + \zr{p,t)\^)] , iV = expf~i / d^p\og{l + \z^{p,t)\^ + \zr{p,t)\^)] . (97) 



Neutrino SU(3) coherent states are characterized by two complex numbers that we denoted by and Zr in Eq. 
(I95|l . We used and z^ for the SU(3) symmetry of the antineutrinos. As a practical convention, we also define 
Ze{p) = Zeip) = 1 in what follows^. 

Evolution operator can be given by the following path integral formula in terms of the SU(3) coherent states: 

{z'{tf),z'{tf)\U\z{U),z{U)) = / (98) 



where the measure is given by Eq. (|127|1 of the Appendix. In this formula 

S[z,z] = J^' dt{zit),zmj^ - Hit)\z{t),zit))-^Hz'(tf),z'itf)\z{tf),z{tf)) (99) 

plays the role of a classical action. The derivation of these formulas and the exact expression for the integral measure 
can be found in the Appendix. 

As before, we write down the Lagrangian^ 

L[z,z]^{^J^-H) (100) 

and solve the Euler-Lagrange equations 

^^^-j^)l[z,z] = and ( ^ A-) L[z,z] ^ (101) 

dt OZn OZn I \dt dZn OZn I 



^ For example 

^|^fe(p)|2 = l + |^^(p)|2 + |2,(p)|2. 
k 

® Here we use the short hand notation (O) = {z{t), z{t)\0\z{t), z{t)), etc. 
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for n = ^.T and those for z* and z*. We use the expectation values of the SU(3) generators which are given below: 

for I, J = e, /i, r (with the convention Ze(p) = ^e(p) = !)■ Here we assumed p € V and p € V. If p ^ V or p ^ V then 
the expectation values are zero. The quadratic terms in the Lagrangian are calculated using the identities 

(T,, ip)T,, {q)) - (T,, (p)) {T,, (q)), {%, {p)T,, (q)) = (p)) (T,, (q)), (103) 

and 

(T,, (p)r,, (q)) = (T,, (p)) (T,, {q)) , (104) 
which are valid for p ^ p and p ^ q. The expectation value of the time derivative term is 



^)=^Y.{1 d'p ^^'^ft + I d'P^T^) (105) 



Using formulas ll()2|l - (|l()5(l we find the Lagrangian as 



Zi(p)z*{p) 



Here Yij [p) is given by 
The equations of motion driven from this Lagrangian are given below: 

iZnip) = X! + ^™ W + Y'mip)) Zi{p) - (7e»(p) + i^ei{p) + Ye.t{p)) Z„(p)zi(p)] , (108) 

i 

iK(P) = E [(7:.(p) - ^n»(p) - y:,(p)) Up) - {iUp) - ^e.(p) - y:m) UpYUp)] , (109) 

i 

where n = fi,T. We write the parameters z^ and Zi in terms of the neutrino wavefunctions as Zj(p) = tpi{p)/ijje{p) and 



and 



.,,fl = fffl = i, = S<|', ->,p, = £ffi. (Ill) 

where the one body wavefunctions are normalized as follows: 

\Mp)\^ + \Mp)\' + \Mp)? = 1, \Mp)\' + \Mp)\' + \Mp)\' = 1- (112) 

If the Eqs. (|110|l and pil|) are substituted in Eqs. p08|l and (|109|l . we find that tpi{p) and ^i(p) obey the following 
Schrodinger equations: 

d ( Mp) \ ( Mp) 

'T, '^^(^) ^{T{p) + ^{p)+Y{p)) i^^p) 1 , (113) 
'^^ \ Mp) / V Mp) 
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d ( MP) \ ( MP) \ 

'Jt ^'^^^^ ^ {-m+m~^y{p)) Mp) ■ (114) 

V ^t(p) / V Mp) I 

Here, the matrix Y is the matrix formed by the elements Yij given in Ea. pU7|) . Note that as one goes from the 
neutrino to antineutrino equations only the signs of the (5m^ terms change. If we substitute Eq. H11U|) and (|lll|l in 
Eq. (|107|l we see that Yij can be written in terms of the one-body wavefunctions as 

= (/ d\R^^i,,{q)M^) - J d^Rp^MQ)rj{q)^ ■ (115) 



VI. CONCLUSIONS 



In this article an algebraic approach to the neutrino propagation in dense media is presented. The Hamiltonian 
describing a gas of neutrinos interacting with each other and background fermions is written in terms of the appropriate 
SU(2) (for two flavors) or SU(3) (for three flavors) operators. Neutrinos as well as antineutrinos are considered. The 
evolution of the resulting many-body problem is formulated as either an SU(2) or an SU(3) coherent-state path integral. 
The evolution operator for the entire system is calculated using two different approximations, namely the saddle-point 
approximation and the operator product linearization approximation. In our case these two approximations yield the 
same answer. This approximate solution of the neutrino evolution is the only one used so far in applications to the 
core-collapse supernovae and the Early Universe. It is important to stress that this solution is only an approximation 
to the many-body problem described by Eq. (|61|l and the corrections to it may play a significant role. 

We should also emphasize that the results for the evolution operator obtained here are applicable to any initial 
state. An initial state such as the one depicted in Eq. Q represents an electron neutrino gas described by a pure 
state, i.e. = p. Such an initial state may be a good description of neutrinos produced in earlier stages of stellar 
evolution, but is inadequate for neutrinos in a thermal distribution. In such cases it would be more appropriate to 
describe the neutrino gas using density operators satisfying the equation 

ip=[H,pl (116) 

where H is the Hamiltonian of Eq. 1)61(1 . The initial density matrix pi, then evolves as p = U piW . Sometimes a 
polarization vector is introduced. In our notation the one-body polarization vector is defined as 

PM) - Tr(J,(g)p) (117) 

for two flavors. For three neutrino flavors one simply calculates the trace above with the SU(3) generators, resulting 
in an eight-dimensional vector. The polarization vector of Eq. ((117|l satisfies the equation 

lP{q)^Tr{[3{q),H]p). (118) 

It is straightforward to show that linearizing the Hamiltonian of Eq. ((61|l using Eq. 144(1 and substituting this 
linearized Hamiltonian in Eq. ((118() yields the evolution equations of the polarization vectors as stated, e.g., in Refs. 
[| and 0. 
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Appendix: Path Integral Representation of the Evolution Operator 

In this appendix, we drive the SU(3) path integral formula in Eq. ((98() . The path integral formula for SU(2) 
coherent states follows very similar lines. 

The coherent states defined in Eq. ((96(1 admit the following resolution of identity: 

[jT 3! d''z^{p)(fz^{p) yr 3! (p) (p) 
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One can use this resolution of identity to write down the path integral formula for the matrix element of the 
evolution operator U. The procedure is well established [s^ and will only be outlined in what follows: We start by 
dividing the time interval [0, T] into N infinitesimally small pieces: 

to^O, ti=e, t2==2e, T^tN^Ne. (120) 

Ignoring the time dependence of the Hamiltonian in each infinitesimal time interval we can write 

(z'(r), z'(T)|C/(r)|z(0), z(0)) = {z'{T), z'(r)|e-*^^(*")e-'^^(*"-i) . . . e^'^^^'^) 1^(0), z(0)). (121) 

We then insert a resolution of identity to the left of each exponential. We will put an additional label to ta to the 
variables Zi{p) and Zi{p) in the resolution of identity which is inserted just next to 6"*"^^^*"^. This way, we obtain 



z'{T),z'{T)\UiT)\z{0),m)^ I TT TT <P^,{p,to.)<fz^{p,t ) ^ ^ 



AT 

X {z\T),z'{T)\z{tN),z{tN)) n (^(^")' z{tc.)\e-''"^'''^Wc.-i),z{to.-i)). (122) 

a=l 

Assuming that only the continuous paths will contribute to the integral when we take N oo limit, we can write 

|z(t„_i), ^(t„-i)) = ( 1 - ) Wo), ziU)). (123) 



In this case the infinitesimal propagator {z{ta), z{ta)\e~^'^^'^*-''^\z{ta~i), z{ta~i)) can be written as 

(z(i„), z(t„)|(l - ieH{t^)){l ~ Sj^Mto.), z{to.)) = e-<^(*-)'^"(*-)l(^^-^(*°»l^(*")'^"(*°». (124) 

Substituting this into Eq. (|122|) and taking the limits N oo and e — > such that Ne = T we arrive at the following 
path integral formula for the propagator: 

(z'(T),z'(T)|C/(r)|z(0),z(0)) = J D[zrzV'^'''^. (125) 

In this formula S[z,z] is given by 

S[z, z]^ dt{z{t), -z{t)\ij^ - H{t)\z{t), -z{t)) - I log(z'(T), Z'(T)|z(r), z(r)) (126) 
and the path integral measure is 

w 3! ^3! 

n n ,^^^y d^ Zf,{p,ta)d^ Zr{p,ta) J| \\_--—n^d^ Z^,{p,ta)d^ Zr{p,ta) ■ (127) 
ot—l p a — 1 p 
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